Anomalous isotopic effect near the charge-ordering quantum criticality 
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Within the Hubbard-Holstein model, we evaluate the various crossover lines marking the opening 
of pseudogaps in the cuprates, which, in our scenario, are ruled by the proximity to a charge-ordering 
quantum criticality (stripe formation). We provide also an analysis of their isotopic dependencies, 
as produced by critical fluctuations. We find no isotopic shift of the temperature T° marked as a 
reduction of the quasiparticle density of states in various experiments, and a substantial positive 
shift of the pseudogap-formation temperature T*. We infer that the superconducting critical tem- 
perature T c has almost no shift in the optimally- and overdoped regimes while it has a small negative 
isotopic shift in the underdoped, which increses upon underdoping. We account also for the possible 
dynamical nature of the charge-ordering transition, and explain in this way the spread of the values 
of T* and of its isotopic shift, obtained with experimental probes having different characteristic 
timescales. 
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There are several experimental evidences |lj| that the 
peculiar properties of the cuprates, both in the normal 
and in the superconducting phase are controlled by a 
Quantum Critical Point (QCP), located near the optimal 
doping 5 = 8 pt- In this framework the phase diagram 
of the cuprates is naturally partitioned into a (nearly) 
ordered, a quantum critical, and a quantum disordered 
region corresponding to the under-, optimally, and over- 
doped regions respectively. The ordered region occurs 
below a second-order transition line T? ($) , which depends 
on the nature of the underlying ordering and, upon in- 
creasing S, ends at T = in a QCP (at 5 — S c > 5 opt )- 
The above correspondence between the theoretical and 
experimental partitioning of the phase diagram, leads to 
a close connection between the hypothetical T?(S) and 
the crossover line T* (S) , below which a pseudogap behav- 
ior is observed in NQR, NMR relaxation-rate, XANES, 
and ARPES measurements [see, e.g., Ref. g] for a re- 
cent overview on La2-aSr,$Cu04 (LSCO)]. In the prox- 
imity of quantum criticality, critical fluctuations medi- 
ate a singular interaction among the quasiparticles which 
can account for both the violation of the Fermi-liquid 
behavior observed in the normal phase of the cuprates, 
and the strong pairing mechanism leading to high-T c su- 
perconductivity [0J|-||] . Various realizations of this sce- 
nario have been proposed, associated with different quan- 
tum criticalities, e.g., antiferromagnetic excitonic [Q, 
change in the symmetry of the superconducting order pa- 
rameter H , or incommensurate charge-density wave || . 

In the QCP framework, any mechanism shifting the 
position of the QCP is mirrored by corresponding shifts 
in T?((5) and in the superconducting critical line T c (5). 
In particular, the observation of isotopic effects (IE's) on 
T c j^U and on T* j|-|ll| suggests that a lattice mech- 
anism underlies the instability marked by the QCP. We 
here consider the single-band Hubbard-Holstein model 
as a minimal model to describe the strongly correlated 
electrons coupled to the lattice, giving rise to a QCP for 



the onset of a phonon-induccd incommensurate charge 
ordering (CO) ||Jl^,|l^]. This introduces density inho- 
mogeneities on a semimicroscopic scale, and corresponds 
to the onset of stripes (m) , coming from the high-doping 
regime. Thus, in our approach, the critical line T?(<5) 
corresponds to the line for CO, Tco(5), which in real 
materials can be masked by pair formation and lattice 
effects. 

In this letter, we first determine the mean-field (m-f) 
critical line T^?q(S), which we identify with the weak- 
pseudogap crossover line T°(8) [3> T*(S)] observed in 
Knight-shift, transport, and static susceptibility mea- 
surements [pj, as the incipient depression of the single- 
particle density of states (DOS). Indeed, for T < T^q 
one expects the CO fluctuations to become substantial, 
leading to a reduction of the quasiparticle DOS, which 
accounts for the weak-pseudogap behavior. Then, we in- 
vestigate the effect of fluctuations near the CO QCP in 
order to determine: a) the fluctuation-corrected critical 
line Tco{&), which we relate to the pseudogap crossover 
line T*(S); b) the IE on T co (S), from which we shall 
also infer the effect on T C (S). In particular we shall de- 
scribe the highly non-trivial effect of quantum criticality 
in determining IE's on T* and T c , which can be strong 
and weak, respectively, and opposite in sign in the un- 
derdoped cuprates. We shall also show why both T* and 
its isotopic shift are observed to be larger in experiments 
with shorter characteristic timescales. The near absence 
of IE on T c near and above 5 op t will also be naturally 
accounted for. 

- The model and the mean-field instability — Our 
two-dimensional (2D) Hamiltonian is 
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where c^J are the fermion operators, a\ are the phonon 
operators, and indicate nearest- and next- 

nearest-neighbor sites, coupled by the hopping parame- 
ters t and t' respectively. The lattice spacing has been set 
to unity. The chemical potential fi is coupled to the local 
electron density = c\ a Ci a , U is the on-site Hubbard 

repulsion, and p q = Ek, CT c k+q,<r Ck ^- ^c(q) (~ Vb/|q| 
at small |q|) is the Coulomb interaction between elec- 
trons in a 2D plane embedded in the three-dimensional 
(3D) space ||,[l3| . loo is the phonon frequency and g is the 
Holstein electron-phonon coupling. In the limit U — > oo, 
this model was solved with a standard slave-boson tech- 
nique, at leading order within a large-N expansion in 
Refs. ||[l3||. The most relevant result was that a CO 
instability with a finite wavevector q c , incommensurate 
with the underlying lattice, was found at T = when 
8 is reduced below a critical value 8c . Besides a small, 
weakly momentum-dependent, residual repulsion, nearby 
this instability the critical charge fluctuations mediate a 
singular scattering 
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of strength V, between the quasiparticles. Here cu n is a 
bosonic Matsubara frequency, 7 ~ t~ l is a characteristic 
timescale, and £q~ 2 is the m-f inverse square correlation 
length, which measures the distance from criticality. Eq. 
(0) displays the behavior of a Gaussian QCP with a dy- 
namical critical index z = 2 fl5| . 

Quite remarkably, the complicated formal structure of 
the quasiparticle scattering, mediated by slave bosons, 
phonons and by the Coulomb interaction, is well repre- 
sented near criticality by a RPA resummation T(q, u> n ) rs 
V e //(q)/[l + V r e //(q)n(q, uj n )\ of an effective static in- 
teraction V e //(q) = U(q) + Vc(q) - Xt @. Here 
n is the fermionic polarization bubble in a 2D lattice, 
A = 2g 2 /tu!o is the dimensionless electron-phonon cou- 
pling, and t/(q) — A + B\q\ 2 is the residual short-range 
repulsion between quasiparticles. For the correspondence 
of V e f f (q) with the parameters of Eq. (P see Ref . [jl6| . 
The instability condition, which occurs for reasonable 
values A - 1, is 1 + V e //(q c )n(q c , u = 0) = 0. At 
T = this determines q c and the position of the m-f 
QCP 8i°\ For realistic parameters we find q c i=s (±1,0) 
or (0,±1), and 8i° ] w 0.2 |,||. By expanding l + V eff U 
near the instability at T = 0, we find £q~ 2 oc 8 — 8c , 
which gives the Gaussian index v = 1/2 for £o- The m-f 
critical line T^q(5) (the solid-line curve in Fig. 1), which 
starts from the m-f QCP at Si°^ is obtained by consid- 
ering the T dependence of the bare polarization bubble, 
which reduces to the simple Fermi- liquid form, oc T 2 , at 



low T. We identify this m-f transition with the exper- 
imental crossover line T°, extrapolating to low T at a 
doping, which we identify with our <5c°^. From the data 
reported in Ref. f§, for LSCO we estimate <^ 0) » 0.22. 
Then we evaluate Tqq(5) by taking standard quasipar- 
ticle (i.e., dressed by the slave bosons) band parameters 
(tgp = 0.2 eV, t' qp = -0.05 eV, u a = 0.07 eV, leading to 
A = 0.2 eV and B = 0.17 eV in the residual repulsion 
U). Vc = 0.22 eV and g — 0.21 eV are adjusted to match 
the experimental extrapolation of T° with the T = m-f 
instability (i.e. 8^)- We point out that the agreement 
between T^l{8) and T°(8) at finite T is obtained with- 
out any further adjustment of the parameters. Similar 
parameters are taken for Bi2Sr2CaCu208 (Bi2212) to fit 
the data in Refs. (T^JH, with V c = 0.22 eV and g = 0.23 
eV. 
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FIG. 1. The phase diagram of the cuprates according to the 
CO-QCP scenario for LSCO (a) and Bi2212 (b). The solid 
line is the m-f critical line ending at T — in the m-f QCP at 
5c ■ The lowest dashed line in panel (a) marks the 3D criti- 
cal line in the presence of fluctuations, ending in the QCP at 
5 C . We took = 1 meV (see text). The dot-dashed line in 
panels (a,b) indicates the "dynamical instability" condition 
(see text) for u) pro ^ e = 1 meV. The intermediate dashed line 
in panel (a) represents the "dynamical instability" condition 
for OJprobe = 1 A* e V- The experimental points for T° (*) and 
for T* measured with fast (x) and slow (□) probes for LSCO 
are from Ref. j|, those for Bi2212 are from Refs. The 
experimental critical temperatures T c are also shown (+)• 

— The phase diagram beyond mean-field — The fluctu- 
ations shift the m-f QCP and critical line to their actual 
position |l!^,^(J. Specifically the fluctuations which me- 
diate the effective critical interaction, Eq. (Q), can be 
included in the polarization bubble via the diagrams of 
Fig. 2a, leading to corrections beyond RPA. From the 
explicit evaluation of these diagrams we find the (self- 
consistent) correction to the mass term m = 7 -1 £~ 2 of 
the fluctuation propagator, Eq. (0), 
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where D = (fi q + |cj„|) 1 is the charge-fluctuation crit- 
ical propagator, u ~ V^ 2 /(7 2 i 3 ) is the coupling resulting 
from the two 4-leg vertices represented in Fig. 2b. The 
critical modes have a 2D dispersion S! q = 7 _1 |q — q c | 2 + 
to, up to an ultraviolet bandwidth cutoff Q m ax ~ 
resulting from the underlying lattice. Since we deal with 
a phonon-driven CO instability, as it also follows from 
the detailed dynamical analysis of Ref. [p"3f , ujq appears 
as the ultraviolet frequency cut-off. By introducing the 
DOS N(Q) = J2q ~ fi q) and tne spectral-density 
representation of D, we rewrite Eq. ([|) in the form 
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where b(z) = [exp(z/T) — 1] 1 is the Bose function. At 
T = 0, Eq. (||) with to = leads to a finite shift of the 
2D QCP S^ 0) -S c (X uq (see Fig. 1). At finite T, as a 
consequence of a constant DOS for the 2D modes, the in- 
tegral in Eq. (0) is logarithmically divergent for to — > 0, 
leading to a vanishing of the renormalized Too (Mermin- 
Wagner theorem). This divergency is removed by consid- 
ering the more realistic anisotropic 3D character of the 
critical fluctuations, introducing a small energy scale 
below which the mode DOS is no longer constant, and 
displays a 3D square-root behavior N(Q < fl±) ~ yfi 
at criticality. This is enough to make the integral in Eq. 
(||) convergent and allows to determine the critical line 
(Eq. ([|) with to = 0) in the anisotropic 3D case as re- 
ported in Fig. la (lowest dashed line). The inclusion of 
fluctuations brings the critical line from temperatures of 
the order of typical electronic energies {T { co ~ t) down 
to much lower temperatures Tco of the order of the ob- 
served T*'s. Indeed, within the CO-QCP scenario, the 
pseudogap arises at T < T* because the quasiparticles 
feel an increasingly strong attractive interaction by ap- 
proaching the critical line Tco{8)- 111 the particle-hole 
channel, this interaction can produce a gap due to the 
incipient CO. At the same time in the particle-particle 
channel, the strong attraction can lead to pair formation 
even in the absence of phase coherence. Therefore T*(5) 
closely tracks the underlying transition line Tco (8). 

On the other hand, the spread in the measured values 
of T* depending on the experimental probe (see, e.g., Ref. 
0), indicates that the CO instability may be "dynami- 
cal", tt± being smaller than (or comparable to) typical 
frequencies of the experimental probes, oj pro b e - In this 
case the self-consistency condition (^) includes uj pro b e as 
the infrared cut-off in the integral over f2, to being re- 
placed by Uprobe- This determines the doping and tem- 
perature dependence of the dynamical instability lines. 



Two examples are reported in Fig. 1 for uj pro b e = lmeV 
(dot-dashed line if Fig. la,b), as in typical neutron scat- 
tering experiments, and uJ W obe — 1 /-teV (second dashed 
line from bottom in Fig. la), as in static experiments 
(NQR, NMR). The corresponding experimental data for 
T* in LSCO are also reported for comparison. We de- 
termine the coupling V between the charge fluctuations 
and the quasiparticles, which is the only parameter for 
which an a priori estimate is difficult, by imposing that 
the fluctuation-corrected QCP is located at the T — ex- 
trapolation of the T*(6) curves. We used 7 = 0.7 eV" 1 
and 7 = 0.4 eV _1 for LSCO and Bi2212 respectively, and 
V = 0.54 eV for both. It is worth noticing that, simi- 
larly to the case of Tqq(S), the agreement between the 
calculated Tco{5) and the experimental points T*{8) is 
obtained without further adjustable parameters. We also 
notice that, contrary to the shift of S c at T = 0, the slope 
of the curve Tco (5) is weakly dependent on ujq. 
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FIG. 2. (a) The two (vertex and selfenergy) corrections 
to the fermionic bubbles (solid line) due to critical charge 
fluctuations (wavy line), (b) Effective 4-leg vertices for the 
critical charge fluctuation field (wavy line) resulting from the 
integration over the fermion loops (solid line). 

- Novel isotopic effects — The m-f weak-pseudogap 
crossover temperature Tqq ~ T° is determined by A only. 
Therefore, it does not depend on uj and is not expected 
to display any isotopic dependence. On the other hand, 
quantities determined by the fluctuations crucially in- 
volve ujq. New physical effects can then arise in the iso- 
topic substitutions. Since ujq decreases by increasing the 



ionic mass, the corrections to 5^ and to T c "^{5) become 
smaller and the fluctuation-corrected quantities 5 C and 
Tco shift to higher doping, closer to the m-f values, as 
shown in Fig. 3, where we report the line TcoWi calcu- 
lated via Eq. (||), with parameters to fit the T* data of 
LSCO (see above), together with its isotopic shift calcu- 
lated for 1g O — * 18 substitution (i.e., for a five percent 
reduction of luq). Correspondingly we expect that the 
portion of the curve T C (S) on the left of T°(6) is rigidly 
translated along the horizontal axis. The rational behind 
this translation is that the whole physics of these ma- 
terials is essentially determined by the proximity to the 
critical line Tco(S) and to the related QCP. On the other 
hand, on the right of the m-f critical line T^^S) ~ T°(8) 
the fluctuations are small and the physical processes at 
T > T°(S) are captured by the m-f description, where 
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A (not luq) is relevant. As a consequence T°(8) and the 
portion of T c {8) near and above T°(5) are not expected 
to be shifted by IE's. 

In the underdoped region, there are two evident con- 
sequences of the isotopic shift, which becomes more sub- 
stantial: The shift upon reducing ujq is negative (as usual) 
in T c , but, contrary to standard theories based on CO 
pseudogap it is positive in Tco ~ T*. Moreover, 
when the slope of Tco{$) is large, a rather small isotopic 
shift in 8 C can result in a substantial shift in Tco ~ T*. 
The steeper is T* , the larger is the IE. On the other hand, 
since the curve T c {8) is rather flat, particularly in the op- 
timal and moderately under-doped regimes, the expected 
IE on T c in these compounds is small in agreement with 
long-standing experiments This large difference in 

the IE for T c and T*, (AT C /AM) / (AT* /AM) < 1, is 
indeed experimentally observed in HoBa2Cu40s (HBCO- 
124), and reported in Ref. |0, where it is also noticed 
that there is a "striking similarity between isotopic sub- 
stitution and underdoping with respect to both T c and 
T*" . Although we are not aware of any systematic analy- 
sis of the doping dependencies of T c , T* and their isotopic 
shifts in HBCO-124, this observation finds its natural in- 
terpretation within our QCP scenario, where the isotopic 
substitution produces a shift of the QCP and is therefore 
nearly equivalent to underdoping. 
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FIG. 3. Calculated effect of the isotopic change le O O 
(i.e. loo = 70 meV and lo' = 66 meV) on T* in LSCO, both 
for a fast- and a slow-probe measurement. The inferred shift 
on T c is also reported 

The slope of Tco{8) increases by increasing Lu pro be, 
while the QCP is unshifted. Therefore, the IE on Tco ~ 
T* is enhanced, and we have the general trend that faster 
probes should detect a larger IE on T*, since the effect of 
fluctuation diminishes. Although we cannot account for 
the near-absent or negative IE on T* within the almost 
static probes in YBa 2 Cu40g (^Jl^l, this general trend 
is in qualitative agreement with a much stronger effect 
observed in the isostructural HBCO-124 with fast neu- 
tron scattering (JH). Indeed, this fast-probe experiment 
should be represented by the curve Tco{5) correspond- 



ing to LJprobe = 1 meV in Fig. 3. The huge IE on T* 
observed in HBCO-124 suggests that the curve Tco{5) 
in this system is steeper than in LSCO. 

Upon underdoping, while T c decreases, its IE incre- 
ses in agreement with a general experimental trend [Q. 
However, in the strongly underdoped materials, where 
a substantial IE in the penetration depth |Q and in 
the x-ray absorption p2] ] have been observed, our re- 
sults can provide only qualitative indications, as other ef- 
fects (magnetic, polaronic, lattice-pinning) not included 
in ([!]) become relevant. For this doping regime, a dif- 
ferent (additional) interpretation of the IE in terms of a 
superconductor-to-insulator QCP has been recently pro- 
posed || . 

Acknowledgments. We acknowledge stimulating dis- 
cussions with C. Castellani and K. A. Muller. 



[1] 

~ 

[4] 
[5] 

[6] 

[7] 



[8] 
[9] 
[10] 

[11] 

[12] 
[13] 
[14] 



[15] 



[16] 
[17] 
[18] 
[19] 



For a summary of this evidences see, e.g., the Introduc- 
tion of C. Castellani, C. Di Castro, and M. Grilli, Z. fiir 



Physik 103, 137 (1997); J. 
Sol. (b) 215, 531 (1999); J. 
Physica C 349, 53 (2001). 
M. Gutmann, E. S. Bozin, 



Tallon, 
Tallon, 



et ai, Phys. Stat, 
and J. W. Loram, 



mat/0009141 



and S. J. L. Billinge, cond- 



A. J. Millis, H. Monien, and D. Pines, Phys. Rev. B 42, 
167 (1990); P. Monthoux, A. V. Balatsky, and D. Pines, 
Phys. Rev. B 46, 14 803 (1992); S. Sachdev and J. Ye, 
Phys. Rev. Lett. 69, 2411 (1992). 
C. M. Varma, Phys. Rev. Lett. 83, 3538 (1999). 

C. Castellani, C. Di Castro, and M. Grilli, Phys. Rev. 
Lett. 75, 4650 (1995). 

M. Vojta, Y. Zhang, and S. Sachdev, Phys. Rev. Lett. 
85, 4940 (2000). 

J. P. Franck, in Physical Properties of High Temperature 
Superconductors IV, edited by D. M. Ginsberg (World 
Scientific, Singapore, 1994), p. 189. 
M. K. Crawford, et ai, Phys. Rev. B 41, 282 (1990). 

F. Raffa, et al, Phys. Rev. Lett. 81, 5912 (1998). 

G. V. M. Williams, et al, Phys. Rev. B 61, R9257 (2000). 

D. Rubio Temprano, et ai, Phys. Rev. Lett. 84, 1990 
(2000). 

A. Perali, et ai, Phys. Rev. B 54, 16 216 (1996). 

F. Becca, et al, Phys. Rev. B 54, 12 443 (1996). 

In our framework the CO instability is a second-order 
transition leading to a gradual increase of the charge 
modulation. It is only when one enters deeply inside the 
(locally) ordered phase that anharmonic distortions of 
the charge profile may arise from the enhanced interac- 
tions with the spin and the lattice degrees of freedom. 
This may lead to the stripe formation, as revealed in 
neutron scattering, see J. M. Tranquada, et al, Nature 
375, 561 (1995). 

This feature is the natural consequence of the damping 
of the charge fluctuations due to the creation of particle- 
hole pairs at arbitrarily low energy when q c joins two 
"hot" points on the Fermi surface. 

G. Seibold, et ai, Eur. Phys. J. B 13, 87 (2000). 

H. Ding, et ai, J. Phys. Chem. Solids 59, 1888 (1998). 
K. Ishida, et al, Phys. Rev. B 58, R5960 (1998). 

J. A. Hertz, Phys. Rev. B 14, 1165 (1976). 



4 



[20] A. J. Millis, Phys. Rev. B 48, 7183 (1993). 

[21] I. Eremin, et al, Phys. Rev. B 56, 11 305 (1997). 

[22] A. Lanzara, et al., J. Phys.: Condens. Matter 11, 



L541 



[23] 
[24] 



(1999). 
J. Hofer, 



et al, Phys. Rev. Lett. 84, 4192 (2000 ) 



T. Schneider and H. Keller, 



:ond-mat/0011381 



5 



